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1 Introduction 

The importance of the mathematical description of correlations in numerous problems of the 
modern statistical mechanics is well-known. Among them in particular, we refer to such fun- 
damental problems as the problem of quantum measurements and of a description of collective 
behavior of interacting particles by quantum kinetic equations pTflO] . Owing to the intrinsic 
complexity and richness of these problems, primarily it is necessary to develop an adequate 
mathematical theory of underlying evolution equations. 

The goal of this paper is to derive rigorously the evolution equations for marginal correlation 
operators that give an equivalent approach to the description of the evolution of states in 
comparison with marginal density operators governed by the quantum BBGKY hierarchy and 
to construct a solution of the corresponding Cauchy problem. 

We briefly outline the results and structure of the paper. In introductory section 2 we set 
forth an approach to the description of the evolution of states of quantum many-particle systems 
within the framework of correlation operators governed by the von Neumann hierarchy [Tl~lll2j . 
In section 3 we introduce the notion of marginal correlation operators. To justify this notion, 
we discuss in detail the motivation of the description of states within the framework of marginal 
correlation operators or in other words, the origin of the microscopic description of correlations 
in quantum many-particle systems. Then we rigorously derive the nonlinear quantum BBGKY 
hierarchy for marginal correlation operators from the von Neumann hierarchy for correlation 
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operators. The nonlinear quantum BBGKY hierarchy gives an alternative method of the de- 
scription of the evolution of states of infinitely many particles in comparison with the quantum 
BBGKY hierarchy for the marginal density operators [13j[Tl]. In section 4 we construct the 
nonperturbative solution of the Cauchy problem of the nonlinear quantum BBGKY hierarchy. 
The nonperturbative solution is determined in the form of expansion over particle clusters, 
which evolution is governed by the corresponding-order cumulant of the nonlinear group of 
operators generated by the solution of the von Neumann hierarchy. The existence theorem for 
initial data from the space of trace-class operators is proved. We also give some comments on 
the mean field (self-consistent field) asymptotic behavior of the constructed solution. Finally, 
in section 5 we conclude with some observations and perspectives for future research in the 
light of the results we present here. 

2 Preliminary facts: the von Neumann hierarchy 

We consider a quantum system of a non-fixed, i.e. arbitrary but finite, number of identical 
(spinless) particles with unit mass m = 1 in the space W,u > 1, that obey the Maxwell- 
Boltzmann statistics. Let = ®^L ^™ be ^ ne Fock space over the Hilbert space 7-L, where 
the n-particle Hilbert space "H n = "H®" is a tensor product of n Hilbert spaces H and we 
adopt the usual convention that "Ho = C. The Hamiltonian H n of the n-particle system is a 
self-adjoint operator with the domain D(H n ) C T-L n 

n n 
1=1 ii<i2 = l 

where K(i) is the operator of a kinetic energy of the % particle and i%) is the operator of 
a two-body interaction potential. In particular on functions ip n that belong to the subspace 
Lq(M. w ) C T>(H n ) C L 2 (M? n ) of infinitely differentiate symmetric functions with compact 
supports the operator K(i) acts according to the the formula: K(i)ip n = —^A qx ip n , where 2nh 
is a Planck constant, and for the operator $ we have: <&(ii,i 2 )'ip n = ^(q^, q i2 )ip n , respectively. 
We assume that the function ^(q^ , is symmetric with respect to permutations of arguments 
and it is translation-invariant bounded function. 

States of a system of the Maxwell-Boltzmann particles belong to the space = 
©^i £ 1 (7/ n ) of sequences / = (/ , fi, ...,/„,. . .) of trace-class operators f n = / n (l, . . . , n) e 
£ 1 ('H n ) and f G C, that satisfy the symmetry condition: f n (l, . . . ,n) = f n (ii, . . . ,i n ) for 
arbitrary . . . , i n ) G (1, . . . , n), equipped with the norm 

oo oo 

II/IUh^h) = ^2\\fn\\zHH n ) = Y1 Tr lv .. jn |/„(l,...,n)|, 

n=0 n=0 

where Tr^...^ are partial traces over 1, . . . , n particles [H]. We denote by ^(J 7 -^) the everywhere 
dense set in 2}(F H ) of finite sequences of degenerate operators with infinitely differentiable 
kernels with compact supports. 

We describe states of a system by means of sequences g(t) = (go, gi(t, 1), . . . , g n (t, 1, . . . , n), 
. . .) G £ 1 (J r -^) of the correlation operators g n {t), n > 1. The evolution of all possible states is 
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determined by the initial-value problem of the von Neumann hierarchy [TT|[T2] 

j t g s {t,Y)=M{Y\g{t)), (1) 

g s (t,Y)\ t=0 = g s (0,Y), s>l, (2) 
where the following notations are used: 

N{Y | g(t)) = -Af s (Y)g s (t,Y) + (3) 

+ Yl J2(-Nint{h,i2))g\x 1 \(t,X 1 )g lX2l {t,X 2 ), 
P:Y=X!\jx 2 heXii 2 eX2 

J2p y=x 1 \jx 2 * s ^ ne sum over a ^ possible partitions P of the set Y = (l,...,s) into two 
nonempty mutually disjoint subsets X\ C Y and X 2 C Y, the operator (— M s ) defined on 
£o(H s ) by the formula 

(-K(Y))f s = ~(H s f s - f.H a ), (4) 
is the generator of the von Neumann equation [15] and the operator (— A/jnt) is defined by 

(-Nlnt(ii,i2))f s = -lmh,i 2 )f s - f s ®(ii,i2))- (5) 
a 

Hereafter we use the following notations: ({Jfi}, . . . , {X|p|}) is a set, elements of which 

are |P| mutually disjoint subsets Xi C Y = (l,...,s) of the partition P : Y = uj^JjXj, 
i.e. |({Xi}, . . . , {X|p|})| = |P|. In view of these notations we state that ({F}) is the set 
consisting of one element Y = (l,...,s) of the partition P (|P| = 1) and |({^})| = 1. We 
introduce the declasterization mapping 8 : ({Xi}, . . . , {X|p|}) —> Y, by the following formula: 
9({Xi}, {X\ P] }) = Y. For example, let X = (1, . . . , s + n), then for the set ({Y}, X\Y) it 
holds: 9{{Y},X\Y) =X. 

On the space £ 1 ('H n ) we also introduce the mapping: M 3 t — > Q n (—t)fn, which is generated 
by the solution of the von Neumann equation of n particles [T5"|[T(3] 

Gn(-t)f n = e-^f n e^. (6) 

This mapping is an isometric strongly continuous group that preserves positivity and self- 
adjointness of operators [T5]. On £,!,(%„) C V(—Af n ) the infinitesimal generator of group §6§ is 
determined by operator (jlj). 

A solution of the Cauchy problem ([I])-© is given by the following expansion [HHT2], [17] 

g s (t,Y) = g(t;Y\g(0)) = (7) 
= %\r\(t,{Xi},.-.,{X\ P \}) J] WO,*), *>1, 

P:Y={JiXi X,CP 

where ^2 P . Y =\j. x x i s ^ ne sum over a ^ possible partitions P of the set Y = (1, . . . , s) into |P| 
nonempty mutually disjoint subsets Xi C Y, the evolution operator 2l|p|(t) is the |P|t/i-order 
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cumulant of groups of operators (J6]) defined by the formula 

2l| P |(t,{X 1 },...,{X|p|})= (8) 
£ (-l) |P ' | - 1 (|P'|-l)! J] G\e(z k )\(-tAZk)). 

P':({Xi},...,{Jf, P ,})=UfcZ* Z k cP' 

Here X)p'-({Xi} {^|p|})=U fc z * * s ^ e sum over a ^ P oss ible partitions P' of the set ({Xl}, . . . , 
{X\p\}) into |P'| nonempty mutually disjoint subsets C ({Xl}, • • • , {X| P |}). For operators 
(J7J) the estimate holds 

\\9s(t)\\ &1{ns) <s\e 2s c% (9) 

where c = maxp : y =UXi (||^ Xl |(0)|| £ i (H|Xi|) , . . . , \\g\x {Pl \(0)\\^(n lX{F{{ ))- 

If <? n (0) G £j("H n ) C -C 1 ('H n ), n > 1, expansion (J7J) is a strong (classical) solution of the 
Cauchy problem ([I])-© and for arbitrary initial data g n (0) G £ 1 (K n ), n > 1, it is a weak 
(generalized) solution [TT|[T2]. 

In case of the absence of correlations between particles at initial time, i.e. initial data 
satisfying a chaos condition, the sequence of correlation operators has the form 

g(0) = (0,^(0, 1),0,...). (10) 

The corresponding solution of the initial- value problem of the von Neumann hierarchy is given 
by the expansion 

s 

g s (t,Y)=Ql s (t,Y)l[g 1 (0,t), (11) 

i=l 

where 2l s (i) is the sth-order cumulant defined by 

% s (t,Y) = £ (-ly^dPI - 1)! II e\*\(-t>Xi)- 

For operators ffTTl) estimate (JH]) takes the corresponding form: 

IW*)L(«.) < s\e s \\ gi (0)\\l 1{Hy 

We note that correlations created in evolutionary process of a system are described by 
formula ffTTl) and determined by the corresponding-order cumulant of the groups of operators 
(J6]) of the von Neumann equations. 



3 The nonlinear quantum BBGKY hierarchy 

The evolution of states of infinite-particle quantum systems is traditionally described by the 
marginal (or s-particle) density operators governed by the quantum BBGKY hierarchy [T51 IT%]. 
In this section we introduce the marginal correlation operators that give an equivalent approach 
to the description of the evolution of such states and describe the nonequilibrium correlations 
in quantum systems. We also rigorously derive the nonlinear quantum BBGKY hierarchy for 
marginal correlation operators from the von Neumann hierarchy (JTJ) for correlation operators. 
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3.1 Marginal correlation operators and marginal density operators 

In the capacity of an example of a mean-value functional of observables [12] we consider the 
definition of the mean-value functional of the additive-type observable = (0, ai(l), . . . , 

Er 1= i a i(^)---) 

oo ^ 

(AW)(t) = V - Tti,...,i+n ai(l)^i+n(*, 1, ■ ■ ■ , 1 + n), (12) 

n=0 

where the operators gi+ n {t), n > 0, are determined by expansions ([7]), and the functional of 
the dispersion for this type of observables 

((A® - (A<»))*)(t) = (13) 
= E n\ Tri '-' 1+ " faiM " (^) 2 {t)) 9l + n {t, 1, . . . , 1 + n) + 

n=0 
oo j 

— r Tri,...,2+n ai(l)ai(2)p 2 +n(*, 1, . . . , 2 + n) . 

n=0 n ' 

For G £(JW) and 5 G £ 1 {F U ) functional (Pjl. fllSjl exists. 

Following to formula (fl3|) . we introduce the marginal correlation operators by the series 

G s (t, 1, . . . , s) = } — Tr a+ i a+n g s + n {t, 1, • • • , s + n), s > 1, (14) 

n=0 

where the sequence g s+n {t, 1, . . . , s + n), n > 0, is a solution of the Cauchy problem of the von 
Neumann hierarchy ([T]). According to estimate 0, series (|14|) exists and the estimate holds: 
||Cr a (t) \\ £ ir n \ < s!(2e 2 ) s c s ^^L (2e 2 ) ri c' 1 . Thus, macroscopic characteristics of fluctuations of 
observables are determined by marginal correlation operators (fT4l on the microscopic level 

({A® - (A«)) 2 )(t) = Tn (a 2 (l) - (^ 1 )) 2 (t))G 1 (t, 1) + Tr li2 o 1 (l)o 1 (2)Gf 2 (t, 1,2). 

Traditionally marginal correlation operators are introduced by means of the cluster expan- 
sions of the marginal density operators F s (t), s > 1, governed by the quantum BBGKY hier- 
archy [13] 

F s (t,Y)= Yi n^w*'**)' ( is ) 

P:Y = \J i X i X,cP 

where Y1p-y=\J Xi * s the sum over a ^ possible partitions P of the set Y = (1, . . . , s) into |P| 
nonempty mutually disjoint subsets Xi C Y. Hereupon solutions of cluster expansions (fT5j) 

G.(t,Y)= (-l) |Phl (|P| - 1)! II s >^ ( 16 ) 

P:Y = \J i X i X,cP 

are interpreted as the operators that describe correlations of many-particle systems. Thus, 
marginal correlation operators (fl6|) are cumulants (semi-invariants) of the marginal density 
operators. 
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As follows from formula ( fl2l) and its generalization [12] the marginal density operators F s (t) 
are defined in terms of the correlation operators of clusters of particles g^ s \t) = (gi +0 (t, {Y}), . . . , 
gi+ n (t, {Y}, s + 1, . . . , s + n), . . .) by the expansion 

F a {t,Y) = J2-Tr 8+1 _ 8+n g 1+n (t,{Y},s + I,..., s + n), s>l, (17) 

n=0 

where the sequence gi +n (t, {Y}, s + 1, . . . , s + n), n > 0, is a solution of the Cauchy problem of 
the von Neumann hierarchy for correlation operators of particle clusters [12J, namely 



g 1+n (t, {Y}, X\Y)= Q{t- {Y}, X \ Y\g(0)) = (18) 
%\ P \{t,{6(X 1 )} t ...,{d(X\ P \)}) J] gwfrXi), s>l,n>0, 

P:({Y},X\Y)=\JiXi X t cP 

where 2l|p|(i) is the |P|£/i-order cumulant defined by formula (jHJ). According to estimate (J9]), 
series (IT7|) exists and the estimate holds: ||-F S (£) || £ w^ \ < e ^ c Y^=o e 3n c n . We note that every 
term of marginal correlation operator expansion (|14|) is determined by the (s + n)-particle 
correlation operator (J7|) as contrasted to marginal density operator expansion (TT7T) which is 
defined by the (1 + n)-particle correlation operator (ITS]) . 

The correlation operators of particle clusters g^ s \t) = (p 1+0 (t, {Y}), . . . ,gi +n (t, {Y}, X\Y), 
. . .) G £ 1 (©^ =0 'H s+n ) can be expressed in terms of correlation operators of particles (J7J) 

g l+ „(t,{Y}.X\Y)= (19) 

e t-DW-'dPi - 1)! n e n 

P:({y},x\y)=Ui^ x iC p P'^p^HJ,-.^ z h cp' 

In particular case n = 0, i.e. the correlation operator of a cluster of \Y\ particles, these relations 
take the form 

g 1+0 ( tl {Y})= J2 H dlx^Xi). 
P:Y={JiXi X t cP 

By the way we observe on that cluster expansions ( 1T5|) follow from definitions (I14p and ( II 7p in 
consequence of relations ( II 9p between correlation operators of particle clusters and correlation 
operators of particles. 

The marginal (s-particle) density operators (TT7]) are determined by the Cauchy problem of 
the quantum BBGKY hierarchy 



^-F s (t, Y) = -M{Y)F,(t, Y) + J2 TW-AUz, * + l))i^+i(t), (20) 



dt 

i£Y 



F s (t)\ t=0 =F s (0), s>l. (21) 

We remind that usually the marginal density operators F s (t), s > 1, are defined by the well- 
known formula of the nonequilibrium grand canonical ensemble [T8l[T9] in terms of the density 
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operators D = (I, D\(t), . . . , D n (t), . . .) governed by the von Neumann equations (the quantum 
Liouville equation) 

F s (t,Y) = (I,D(t))- 1 ^^Tr a+li ..., a+n D B+n {t,X), 

n=0 

where (I,D(t)) = Y^=on[ r ^ ri ,---,nDn('t) is a normalizing factor, I is the identity operator and 
Y = (1, . . . , s), X = (1, . . . , s + n). Thus, along with the definition within the framework of 
the non-equilibrium grand canonical ensemble the marginal density operators can be defined 
within the framework of dynamics of correlations that allows to give the rigorous meaning of 
the states for more general classes of operators than the trace-class operators. 

If -F(O) G 2} a (Fu) = ©^Lo"" ^ai^n) and a > e, then for t G R a unique solution of 
the Cauchy problem (|20I) - (|2T1) of the quantum BBGKY hierarchy exists and is given by the 
expansion [IT] . [20J 

F s (t, Y) = V - Tr s+1 ,..., s+n 2l 1+ „(t, {Y}, X\Y)F s+n (0, X), s > 1, (22) 

n=0 

where the (1 + n)th-order cumulant 2ti +ri (t) of groups of operators (EJ) is defined by 

% 1+n (t,{Y},X\Y)= Yl (-l) |Phl (|P|-l)! II felH-W). ( 23 ) 

p : ({y}^\^)=U^ ^cp 

J2p is the sum over an possible partitions P of the set {{Y}, X\Y) into |P| nonempty mutually 
disjoint subsets X { C ({Y}, X \ Y). 

Formally, the evolution equations for marginal correlation operators are derived from the 
quantum BBGKY hierarchy for marginal density operators (|20l) on basis of expression ( fl6l) . 
Then the evolution of all possible states of quantum many-particle systems obeying the Maxwell- 
Boltzmann statistics with the Hamiltonian ([TJ can be described within the framework of 
marginal correlation operators governed by the nonlinear quantum BBGKY hierarchy 

jG s {t, Y) = Af(Y | G{t)) + Tr s+1 ^(-JVteft s + 1)) {G s+1 (t, Y,s + 1)+ (24) 

+ Yl G lXll (t,X 1 )G lX2l (t,X 2 )), 

P : (Y,s + 1) =X 1 (JX 2 , 
ie X 1 ;s + le x 2 

G s (t,Y)\ t=Q = G s (0,Y), 8 >1. (25) 

In equation ([24"]) the operator M{Y \ G(t)) is generator of the von Neumann hierarchy (flj) 
defined by formula ([3]), i.e. 

Af{Y\G{t)) = {-M s (Y))G s {t,Y) + 
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where the operators (— AQ and (— A/j n t) are defined by (J4j) and ([5]) respectively, J2p-y=x 1 {Jx 2 
is the sum over all possible partitions P of the set Y = (1, . . . , s) into two nonempty mutually 
disjoint subsets X\ C Y and A 2 C Y, and p-.(y,s + i) = X!\jx 2 , * s the sum over all possible 

i 6 Xi;s + l e x 2 

partitions of the set (Y, s + 1) into two mutually disjoint subsets X\ and X 2 such that ith 
particle belongs to the subset X\ and (s + l)th particle belongs to X 2 . As far as we know 
hierarchy fliH|) was introduced by M.M. Bogolyubov [13] and in the papers of J. Yvon [21J and 
M.S. Green [22] for systems of classical particles. 

Another method of the justification of evolution equations for marginal correlation operators 
consists in their derivation from the von Neumann hierarchy for correlation operators ([I]) on 
the basis of definition (HU). 



3.2 The derivation of the nonlinear quantum BBGKY hierarchy 

In this section we establish that marginal correlation operators (JHJ) are governed by the non- 
linear quantum BBGKY hierarchy With this aim we differentiate by time variable the 
marginal correlation operators defined by series (TH]) in the sense of the pointwise convergence 
on the space £ 1 ('H S ) 

jGsitX) = Xl^TWw+n {(-K+n(X))g s+n (t,X) + 

n=0 

+ E E ^(-AW^2))W^i)W^)), (26) 

P:X=Xiux 2 iieii j 2 ex 2 

where we use the notations: A = (1, . . . , s + n), Y = (1, . . . , s), ^p y=jfiU^2 * s ^ ne sum over 
all possible partitions P of the set Y = (1, . . . , s) into two nonempty mutually disjoint subsets 
Xi C Y and A 2 C Y, and operators (— M s +n) an d (— A/j n t) are defined by formulas (jl]) and §5§ 
respectively. Taking into account the equality 

-M s+n (X) = -M s (Y)-M n (X\Y) + J2 E (-Mnt(ii^ 2 )), 

h£Y i 2 £X\Y 

and the identity 

Tr s+1 ,... iS+ „(-A4(X \ Y))g s+n (t) = 0, 

and according to the symmetry property of the correlation operators g s+n (t, X), for the first 
term in the right-hand side of identity ( 12"6l) in terms of definition ( fl4l) we obtain 



OO j 

^ r Tr s+1) ... iS+n (-AC+n(A))c/ s+n (*, 1, . . . , s + n) = (27) 

n=0 

oo 

-A4(Y)G s (t,Y)+J]-Tr s+w+n J2 (-Mnt(ii,i2))gs+n(t,X) = 
n=i ' heYi 2 eX\Y 

-K(Y)G a (t, Y) + J2 — Tr s +i,..., s+ „+i ^(-A/" int (ii, s + l))c/ s+n+1 (t, X, s + n + 1) = 



n=0 uev 
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-Af s (Y)G s (t, Y) + Tr s+1 ("-Mnt^i, 8 + l))G s+1 (t, 1, . . . , s + 1). 



Let us consider successively the following four parts of the second term of the right-hand side 
of identity ([26} 



E E E<-^..^»= E ( E E + < 28 > 

P:X=X 1 UX 2 heXi i 2 €X 2 P:X=X 1 UX 2 iiGYnXi i 2 GYnX 2 

+ E E + E E + E E )hvu(*i,«). 

nex\ynXi j 2 eynx 2 ueYnXi i 2 ex\Ynx 2 uex\YnXi i 2 ex\Ynx 2 
Taking into account that the equality is true 

E g\x 1 \(t,X 1 )g lX2l (t,X 2 ) = 

P : X = Xi UX 2 , 

E E ^( t . y i.%i+™(*. ir >.(^\ y )\4 

P:Y = Y 1 |JY 2 , ZcX\Y 
y ^ 0, y 2 ^ 

and the validity of the following equality (according to the symmetry property of operators 

9n(t)) 

Tr s+lt ___ jS+n Y 9\Y 1 \+\z\(t,Y 1 ,Z)g\ Y2]+KX \ Y) \z\(t,Y 2 ,(X\Y) \Z) = 

ZCX\Y 

n . 
v — v TV 

= Tr s+lj ... jS+n — - — g lYl \+n-k(t, Y u s + l,...,s + n-k)x 
x g\Y 2 \+k(t, Y 2 , s + n - k + 1, . . . , s + n), 
for the first part of equality (12"5|) of the second term of identity (1251) it holds 

1i , a +l,...,«+n E E E (-■Mnt( i l>«2))^il(*»- X 'l)^Jfa|(*>- X '2) = 

P:X=XiUX 2 nernXi i 2 eYnx 2 



n . 

nl 

x 



= Tr s+1 ,..., s+n ]T E E (-MUh,i2)) E - A ,, fn _ /,), 

P:Y=Y 1 \JY 2 heY 1 i 2 &2 k=0 ' V >' 

xg\Y!\+n-k{t, Y u s + 1, . . . , s + n - k)g\Y 2 \ +k {t, Y 2 , s + n - k + 1, . . . , s + n). 
Then in terms of definition (TH|) the last expression takes the form 



nl 

x 



E^t' 1 ' 1 '- 1 E E E(-- A/ i^ i ' i2 ))E /C ! (// _/, ) ! 

n=o p : y=y 1 uy 2 h&i % 2 & 2 k=o v ; 

xg\ Yl \ +n -k(t, Yi, s + 1, . . . , s + n - k)g\Y 2 \+h(t, Y 2 , s + n - k + 1, . . . , s + n) 

= E E ("A^nt(zi, i 2 ))C|^| (t, l^Cl^lCt, K 2 ). 

P:Y=Y 1 \JY 2 heYii 2 EY 2 
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Hence for the first part of equality (1251) of the second term of the right-hand side of identity 
( |26i) we have 

X 

^rjTWw+n £ E E (-M nt (2 1 ,2 2 ))^ 1 |(t,X 1 )^ 2 |(t,X 2 )= (29) 

n=0 ' P:X=XiUX 2 h&YnXt i 2 &V\X 2 

= E E E^^^'^))^!^^)^!^^). 

Similarly the second and third parts of equality (12"51) of the second term of the right-hand side 
of identity (126]) are expressed in terms of definition (Tl4j) in the form 

oo 

£^+. £ ( £ £ + (») 

n=0 P:X=XiUX 2 iieX\ynXi i 2 eYnX 2 

+ E E )(-' / ^*(*l»*2))^ x |(t,Xl)»[X a |(t,^2) = 

jiernXi j 2 ex\ynx 2 

= ^Tr s+1 (-M nt (z, S + l)) ^ G^i^XOG^i^Xa), 

ier P : (Y,s + 1) = X 1 (JX 2 , 

i g Xi ; s + 1 £ X 2 

where X] p : (y,s + i) = XiMx 2 , * s the sum over a ^ possible partitions of the set (Y,s+ 1) into 

% £ Xi ;s + l£ X 2 

two mutually disjoint subsets X\ and X 2 such that ith particle belongs to the subset Xi and 
(s + l)th particle belongs to X 2 . 

Then taking into account the following identity for the fourth part ( |28|) of the second term 
of ([26]) 

— TT s+i,...,s+n E E E (-M a t(h,i2))g\x 1 \(t,X 1 )g lX2l (t,X 2 ) = 0, 
n=o ' P:X=Xiux 2 ijexvynXi « 2 GX\ynx 2 

and identities (T2"9"1) . (15U|) . for the second term of the right-hand side of (1251) it holds 

E ~j E E E (-• A/ int( i 1^2))5|X 1 |(*,^l)5|X 2 |(i,^ 2 ) = (31) 

n=0 ' P:X=XiUX 2 iiSXi j 2 £X 2 

= E E E(- A/ "^(^i^ 2 ))G| yi |(t,y 1 )G'| y2 |(t,y 2 ) + 

P: Y=Y! |J Y" 2 iiSYi i 2 eF 2 

+ ^Tr s+1 (-jV int (i,s + l)) £ G ]Xll (t,X 1 )G ]X2] (t,X 2 ). 

ieY p : (y,s + i) = Xi|jx 2 , 

i 6 Xi ; s + 1 6 X 2 

In consequence of identities ( 1271) and (l3Tj) we finally derive 

^G a (t,Y) = -M(y)G.(t,Y)+ Yl E E(-^(« 1 ' i2 ) G l y il(*' yi )G|y 3 |(*,y2)+ 

p : y=y x U y 2 iieyi j 2 ey 2 
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+ J2^s+i(-^t(i,s + l))(G s+l (t)+ G lXll {t,X 1 )G lX2l (t,X 2 )), 

i£ Y P : (Y,s + 1) =Xi|J X 2 , 

i e X 1 ; s + 1 G X 2 

where we use notations accepted above in (1261) . (1301) . The constructed identity for the marginal 
correlation operators defined by expansion (|14p we treat as the hierarchy of evolution equations, 
which governs the marginal correlation operators of quantum many-particle systems. 

We also formulate the nonlinear quantum BBGKY hierarchy in case of many-particle systems 
obeying quantum statistics with the Hamiltonian 

n n n 

H n = Y,m + Y. E * w («i,-...«i*), 

i=l k=l ii<...<ifc=l 

where $^ is a fc-body interaction potential. We introduce the symmetrization operator 
and the anti-symmetrization operator S~ on "H®" which are defined on the space 1-L n by the 
formula 



5 « =^E( ±1 ) M ^ 



where the operator p n is a transposition operator of the permutation 7r from the permutation 
group & n of the set (1, . . . , n) and \tt\ denotes the number of transpositions in the permutation 
7r. The operators S„ are orthogonal projectors, i.e. (S^) 2 = S^, ranges of which are corre- 
spondingly the symmetric tensor product and the antisymmetric tensor product H~ of n 
Hilbert spaces Tl. We denote by J 7 ^ = ©^Lq'H^ an d = ®^=o^-n the Bose and Fermi Fock 
spaces over the Hilbert space % respectively. 

The evolution of all possible states of many-particle systems of bosons or fermions is described 
within the framework of marginal correlation operators G(t) = (Go, Gi(t), G s (t), .. .) G 
£ 1 (J 7 ^) governed by the following nonlinear quantum BBGKY hierarchy 

jG s {t, Y) = -M s (Y)G s (t, Y) + (32) 

p 

+ E E ••• E {-MF Zr \zx,...,z lPl ))sti[G lXil (t,x i ) + 

P:Y = [J i X i , Z x Cli, Z|p|CX| P |, X S CP 

|P| > 1 2i^l Z| P | ^ 

oo n s 

+ EE E Tr s+li ... iS+1+n _ fe (-^ +1) (ji,...,j fc ,s + l,...,s + l + n- k)) x 

n=l fc=l l=ji<...<i fc 

e ^w* n <wf,*,), 

P:(l,...,s + l + n-fc) = U=i*i> ^iCP 
^P\(ji,-,JU |P|<n + l 

where notations accepted above are used and the operator (—Af^) acts on £q(H s ) C £ l (T-L s ) 
according to the formula 

(-■A/ffft, • • • , i k ))f s = -^($ (fc) (*i, • • • , i k )f s - fs& k) ^i, ■ ■ • , t*)). (33) 
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We emphasize that the evolution of marginal correlation operators of both finitely and in- 
finitely many quantum particles is described by initial- value problem of the nonlinear quantum 
BBGKY hierarchy (|24p (or (1321) ). For finitely many particles the nonlinear quantum BBGKY 
hierarchy is equivalent to the von Neumann hierarchy ([I]). 

4 A nonperturbative solution of the nonlinear quantum 
BBGKY hierarchy 

In this section we construct a solution of the initial-value problem of the nonlinear quantum 
BBGKY hierarchy (1241) . A nonperturbative solution is represented in form of an expansion 
over particle clusters which evolution is governed by the corresponding-order cumulant (semi- 
invariant) of nonlinear groups of operators ([7]) generated by the von Neumann hierarchy ([I]). 
A solution representation in the form of the perturbation (iteration) series of hierarchy ( 12~4"|) is 
derived as a result of applying of analogs of the Duhamel equation to cumulants of groups of 
operators ([7]) of constructed solution. 

4.1 A solution in case of chaos initial data 

To construct a nonperturbative solution of the Cauchy problem (1241) - (125]) of the nonlinear 
quantum BBGKY hierarchy we first consider its structure for physically motivated example of 
initial data, namely, initial data satisfying a chaos property 

G s (t,Y)\ t=0 = G 1 (0,l)5 s , u s>l, (34) 

where 5 s .i is a Kronecker symbol. Chaos property (I34p means the absence of state correlations 
in a system at the initial time. 

According to definition (Ti~4"l) and solution expansion (TIT]) , in the case under consideration the 
following relation between the marginal correlation operators and correlation operators is true 

Gxfai) = 9l (0,i). (35) 

Taking into account the form fTTTj) of a solution of the initial-value problem of the von Neumann 
hierarchy ([1]) in case of initial data f fTOj) . for expansion (fT4"j) we obtain 

oo j s+n 

G s (t,Y) = ^-Tr fl+w+n a s+n (*,l,...,s + n) Y[gi(0,i), ( 36 ) 

n=0 ' i=l 

where Vls+nit) is (s + n)t/i-order cumulant ([8]). In consequence of relation (1351) we finally derive 

oo ^ s+n 

G s (t,Y) = ^-jTr s+lr .. )S+n a s+n (t,l,...,s + n) JjGi(0,i), s > 1. (37) 

n=0 ' i=l 

If ||Gi(0)||£i(^) < (2e) _1 , series ( 1371) converges, since for cumulants (IE]) the estimate holds [H] 

l|a»(*)/ll^(w.) < n\e» im| £l( «„). 
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From the structure of series ( |37|) it is clear that in case of absence of correlations at initial instant 
in a system the correlations generated by the dynamics of quantum many-particle systems are 
completely governed by cumulants (|8]) of groups of operators 

Thus, the cumulant structure of solution ([7]) of the von Neumann hierarchy ([1]) induces 
the cumulant structure of solution expansion (|37|) of the initial-value problem of the quantum 
nonlinear BBGKY hierarchy for marginal correlation operators. 

The evolution equations which satisfy expression (l3"Tj) are derived similarly to the derivation 
of hierarchy ff24"j) given in section 3.2 on the base of definition fl36|) . 

We note, that in case of initial data f lTUj) solution ffTTj) of the Cauchy problem (DO)-© of the 
von Neumann hierarchy may be rewritten in another representation. For n — 1, we have 

fll (t,i) = 5ai(t,i)fli(o,i). 

Then, within the context of the definition of the first-order cumulant, 2ti(— t), and the dual 
group of operators 2li(£), we express the correlation operators g s (t), s > 2, in terms of the 
one-particle correlation operator g\{t) using formula f fTTj) . Hence for s > 2 formula ffTTj) is 
represented in the form of the functional with respect to one-particle correlation operators 

s 

g s (t, Y | 9l (t)) = Y) J] 9l {t, i), s>2, 

i=i 

where QL s (t,Y) is sth-order cumulant (JS} of the scattering operators 

s 

Q s (t, Y) = g a (-t, Y) fj i), s > 1. (38) 
i=i 

The generator of the scattering operator Gt(Y) is determined by the operator 

niw = E E H^fffr. • • • > <*)). 

fc=2 H<...<ifc = l 

where the operator (— A/"^) acts on £q("H s ) C £ 1 ('H s ) according to formula (|3"3"|) . 

Similar representation of a solution holds for marginal correlation operators (137)) which forms 
inherently the basis of the kinetic description of the evolution. In this case the marginal 
correlation functional G s (t,Y \ G\(t)), s > 2, are represented by the expansions 

oo . s+n 

G s (t, Y | G 1 (t)) = E ~ Tr s+1) ..., s+n 9J 1+n (£, 0({Y}), s + 1, . . . , s + n) J[ G^t, i), (39) 

n=0 ' i=l 

where the operator Gi(t, i) is given by f[3^|) for s — 1, and we use the notion of the declasteriza- 
tion mapping defined in section 2 [20]. In expansion ( )39|) the (1 + n)t/i-order evolution operator 
%Ji+n{t) is defined by the formula [9] 
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n n n— ni — n fc _! 

awMan,*^)^ £(-!)*£... £ („-„,-...-„»), * 

fc=0 ni=l n fc =l v 1 K; 

x2l s+n _ ni _..._ nfc (t, {Y}, s + l,...,s + n-n 1 -...-n k ) x 

s+n— 111- flj 

X II E loqr E II Tjrn 2ti + ix,,(Mz,,^), 

i=i i>, U V- . 1 ii/...#i| Dj i=i x^cDj 1 y' - 

|Dj I < s + n — ni — ■■• — rij 

where X/Dj-^HJj.-Xt * s the sum over a ^ possible dissections of the linearly ordered set 
Zj = (s + n — ni — ... — nj + l,...,s + n — n 1 — ... — rij-i) on no more than s + n — ri\ — . . . — rij 
linearly ordered subsets, and the operator 2li +n (£) is the (1 + n)-order cumulant ([8]) of the 
scattering operators (1581) . For example, the lower orders evolution operators 5Ji +n (t, 9({Y}), s+ 
1, . . . , s + rij , n > 0, have the form 

*J 1 (t,e({Y})) = % 8 (t,e({Y}), 

s 

gj 2 (t, 0({y», a + 1) = 2t s+1 (t, 0({y», s + 1) - 2t s (t, 0({y») ^ 2t 2 (t, z, a + 1), 

i=l 

and in case of s — 2, it holds 

9J 1 (t,0({l,2})) = &(t,l,2)-J. 

We point out also that in case of chaos initial data solution expansion ([22"]) of the quantum 
BBGKY hierarchy (120]) for marginal density operators differs from solution expansion (|3Tj) of 
the nonlinear quantum BBGKY hierarchy (1241) for marginal correlation operators only by the 
order of the cumulants of the groups of operators of the von Neumann equations [12] , [2D] 

x s+n 

F s (t,Y) = J2^s + i,...,s +n ^i + n(t,{Y},X\Y) J]>i(M, ( 4 °) 

n=0 ' i=l 

where 2li +ri (t) is the (1 + n)t/i-order cumulant fl23|) . Series (j4"Ul) converges under the condition: 
11^(0)11^^) < e- 1 . 



4.2 The structure of a nonperturbative solution expansion 

The direct method of the construction of a solution of the nonlinear quantum BBGKY hi- 
erarchy ( |24"|) in the form of nonperturbative expansion consists in its derivation on the basis 
of expansions ( TT6|) from nonperturbative solution (122"]) of initial-value problem of the quan- 
tum BBGKY hierarchy (12"0]) -( 12"T]) . Following stated above approach, we derive a formula for 
a solution of the quantum nonlinear BBGKY hierarchy for marginal correlation operators in 
case of general initial data on the basis of definition ( fl4l and nonperturbative solution (JTj) of 
initial- value problem of the von Neumann hierarchy (JI])-(J2J)- With this aim on f n 6 2}{1-L n ) we 
introduce an analogue of the annihilation operator 



(o/) s (l, . . . , s) = Tr s+1 / s+1 (l, . . . , s, s + 1), s > 1, 



(41) 
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and, therefore we have 

°° (±l) n 

(e ±a /)s(l, • • • , s ) = / ] j— T^B+l,...,a+nf a +n{h ■ ■ ■ , S + n) . 



n=0 



According to definition (fl4j) of the marginal correlation operators, i.e. 

G(t) = e a g(t), 

where the sequence g(t) is a solution of the von Neumann hierarchy for correlation operators 
defined by group (171), i.e. g(t) = Q{t \ g(0)), and to the equality: g(0) = e~ a G(0), we finally 
derive 

G(t) = e a Q{t I e~ a G(0)). (42) 

To set down formula (T4"2"j) in componentwise form we observe, that the following equality 
holds 

J] (e-'G(0)),*|(*) = ^Tr wl ,..., s+n+fc £ kUk L kV ■ ■ ■ ( 43 ) 

X,cP fc=0 ' fci=0 X >' 

fc |P|-2 , , 

... - — |P| ~ 2 " -G\ Xl \+k-ki{Q,Xi,s + n + l,...,s + n + k- kx) ... 

*-^> fc|p|_i!(fc|p|_ 2 - fc|P|_i)! 

fc| P |_i=0 11 1 ii' 

• • • G |^|p||+A;|p|-i( ' X \ p \> s + n + k ~ k \P\-i + 1, . . . , a + n + fc). 
Then according to formulas fj42|) and ([7]), for s > 1 we have 

00 

G.(t, F) = - { Tr s+ i,..., s+ „. £ 2l| P | (t, {X,}, . . . , {X| P |}) J] ( e - o G(0))| Xi |(X,), 

n=0 ' P:(l,...,s+n)=Ui^ X * cP 

where 2l|p|(t) is |P|£/i-order cumulant ([8]), and as a result of the validity of equality fj43]) for 
sequence (l4"2l we obtain 

G.(t,l,...,s) = 

= E ^ Ti w,..., s+ . E(-i) %,, w %, E %i (*> w» • • • > x 

n=0 fc=0 ^ ' P: (l,...,s+n-fe)=Uj Xi 

k ^, fc |P|-2 ^ , 

s + n-k + l,...,s + n-ki)... G|Xn,||+fc, P |_ x (0, X\ P \,s + n- k\ P \-i + 1, . . . , s + n). 

Consequently the solution expansion of the nonlinear quantum BBGKY hierarchy has the 
following structure 

00 

G s {t,Y) = V-rTr s+w+n [/ 1+n (t;{F}, S + l,... )S + n| G(0)), s > 1, (44) 
^— ' n! 

n=0 
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where we introduce the notion of the (1 + n)th-order reduced cumulant U\ +n (t) of nonlinear 
groups of operators (J7j) 

U 1+n (t; {Y}, s + l,...,s + n | G(0)) = (45) 

fe=0 V ; P:(e({l,...,s}),s+l,...,s+n-k)=\J i X i 

fc » j fc |P|-2 £ J 

x 2 / 7~rn — rrr • • • E i ui 7 ^G\xi\+k-ki(Q,x x , 

z_, _ ^ |P| _ 1 !(fc |P| _ 2 - fcpi-i)! 

s + n- fc + l,...,s + n-A;i)... G|X|p||+fc| P |_ 1 (0 ) X| P |, s + n - fc| P |_i + l,...,s + n). 
We give simplest examples of reduced nonlinear cumulants ( 1451) : 



^(t;{y}|G(o)) = ^ ; r|G(o)) = 

£ ^(^{^....{Xip,}) n G| Xi |(0,X,), 

P:(l,..., a )=U 1 X i XiCP 

[/ 2 (t;{y}, S + l|G(0))= ^ SIipi^^},...,^,}) IJG^CO,^)- 

P:(y,s+l)=Ui^« x » cP 

IPI 

£ 0|P|(t,{X 1 },...,{X| P |})J] J] C7|x 4 |(0 > X 4 )G|x J | + i(0,X, >fl + l). 

P:(l,...,«)=U i ^ i =1 XiCP, 

Xi ^ Xj 



We remark that in case of solution expansion ([22]) of the quantum BBGKY hierarchy, an 
analog of reduced cumulant fj45|) is the reduced cumulant of groups of operators ([6]) defined by 
formula [H] 

*7 1+n (t; {y}, s + 1, . . . , s + n) = ^(-^ k k] ^i k y Gs+n-k(-t). 



k=0 



4.3 Reduced cumulants of nonlinear groups of operators 

We indicate some properties of reduced nonlinear cumulants (j4"5l) of groups of operators ([7]). 
According to formula f|44|) and properties of cumulants (jSJ), namely 2l„(0) = IS n; x, the following 
equality holds 



U 1+n (0;{Y},s + l,...,s + n | G(0)) 
. n! 
fc!(ra- A;)! 



Sti(0, {1, . . . , 8 + n - k})G s+n {0, l,...,s + v) 

k=0 



= G s+n (0, l,...,s + ri)5 nfi , 
and hence the marginal correlation operators determined by series ( l4*4l) satisfy initial data ( 1251) . 
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In case of n = for / e ^(J 7 ^) in the sense of the norm convergence of the space 2}{1-L S ) 
the infinitesimal generator of first-order reduced cumulant ( )45j) coincides with generator ([3]) of 
the von Neumann hierarchy (pQ) 

hm^fM*; {Y} I /) - f s (Y)) = M{Y | /), s > 1, 

where the operator N(Y \ f) is defined by formula ([31). In case of n = 1 for second-order 
reduced cumulant f j4"5l) in the same sense we obtain the following equality 

Tr s+1 lim \ U 2 (t; {Y}, s + 1 | /) = V Ti s+1 (-M mt (i, s + 1)) (f s+1 (t, Y,s + 1) + 



+ f\x 1 \(t,X 1 )f [Xa \(t,X 2 )), 



P:(y,s + l) = Xi|J^2, 
!£Xi;s + leX 2 

where notations are used as above for hierarchy (124"|) . and for n > 2 as a consequence of the 
fact that we consider a system of particles interacting by a two-body potential, it holds 

Tr s+1 s+ „ lim - U 1+n (t; {Y}, s + 1, . . . , s + n | /) = 0. 

In case of initial data satisfying a chaos property, i.e. G^(0) = (0, Gi(0, 1),0, . . .), for the 
(1 + njth-order reduced cumulant we have 

s+n 

U 1+n (t;{Y},s + l,...,s + n \ G {1) (0)) = 2l s+n (t, 1, . . . , s + n) J] Gi(0, i), 

i=i 

i.e. the only summand that gives contribution to the result is the one with k = and |P| = s+n, 
since otherwise there is at least one operator G s (0) with s > 2 in the last product. We note 
that in section 4.1 the same result was obtained using the properties of solution of the von 
Neumann hierarchy ((T]). 

For the (1 + n)th-order reduced cumulant fj45|) the following inequality holds 

\\U 1+n (t; {F}, s+l,...,s + n\ /)|| £l(w<+fJ < 2n\s\(2e 3 y +n c s+n , (46) 



wher6C = P n T aX n iv uu ^ W\X 1 \+k-k l WH lXll+k . hl ),--- 

P: (1, . . . , s + n — k) = \J. Xi k, fci, . . . , fcipi_i £ 

6 (s + n — fc + 1, . . . , s + n) 

H/l^iPil+fciPi-ilUHWix^i+fep^^)- 
To prove this inequality we first remark that for cumulant (jSJ) the following estimate holds 

\\Ql lPl (t,{X 1 },...,{X lPl })f n \\ £l{Hn) < |P|!el p l||/ n || £l(Wn) . (47) 

Indeed, we have 

l|2l|P|(^ {^l}) ■ ■ ■ ) {^|P|}) fn\\ - 

< E (IP'I- 1 )'!! II fe)|(-^^))/n|UH«») = 

p' : ({x 1 },...,{X| P |})=u fe ^ ^ fc cp' 

|p| 

= ||/„|| £l( ^)Es(|PM)(Z-l)!, 
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where s(|P|, I) are the Stirling numbers of second kind and we use the isometric property of the 
groups Q n (—t), n>l. Estimate ( H71) holds as a consequence of the inequality 

|P| 

5>(|P|,Z)(/-l)!<|P|!el p l 
1=1 

Then owing to estimate (j4"T|) . for the (1 + n)t/i-order reduced cumulant (|45p we have 
\\U 1+n (t; {Y},s + l,...,8 + n\ f)\\ zHHa+n) < 

V |P|'e |P| V kl 



k=0 v > P: (l,...,s+n-fc)=Ui X t k 1= v 17 



fc|P|-i=0 



*-^(n-k)\ ^ 

k=0 V ' P:(l,...,s+n-k)=[j l X l 

As result of using of the definition of the Stirling numbers of second kind s(s + n — k, I) and 
the inequalities 

n . n . s+n— fc 

n! oiDi i n! 



E^ E ipn 2|p| -'-E(^ E»(« + »-*.')«e»-< 

fc=0 V 7 P: (l,...,s+n-Jfc)=Ui X » fc=0 /=1 

^ d( ' + _ n ~, fc)! e 3(s+ "' fc) < 2n!,!(2e 3 ) s+ ", 



n 



we obtain estimate ( 146]) . 

Thus, according to estimate f|4"6j) . for initial data from the space £ 1 ('H n ) series (jH]) converges 
under the condition: c = max n >i ||G n (0) || . < (2c 3 )" 1 , and the following inequality holds 



n=0 



A solution of the Cauchy problem of the nonlinear quantum BBGKY hierarchy for marginal 
correlation operators (l24"j) is determined by the following one-parametric mapping 



R 3 t -> Uit | /) = e n £(t | e -0 /), (49) 

which is defined on the space £ 1 (J 7 ^) owing to estimate (JlSj) . and has the group property 

U{h\U{t 2 | /)) = U(t 2 \U{h I f))=U(h + t 2 i /). 

Indeed, according to definition (14"T|) and taking into consideration the group property of the 
mapping Q{t | ■), we obtain 

U(h + t 2 | /) = e a 0(t! + t 2 | e"V) = e*g{t x \ Q{t 2 \ e~*f) = 

= e a G{h | e~ a e a g(t 2 \ e~ a f) = e a Q{t x \ e~ a U(t 2 | /)) = U(t x \ U{t 2 | /)). 
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To construct the generator of the strong continuous group U(t;Y | •) we differentiate it in 
the sense of the norm convergence on the space £ 1 ('H S ) 

jUfrYl f)\t=o = j t (e a G(t | e~ a f)) s (Y)\ t=0 = 

oo 1 

= E~i Tr s+ W +n JV(X I G(t I e-°/))| t =o = (e«W(- I e~ a f)) s (Y), 

n=0 

where jV(- | /) is a generator of the von Neumann hierarchy ([1]) defined by formula ([3]) on the 
subspaces £q("H s ) C -C 1 ('H s ), s > 1, or in the componentwise form 

(e a Af(- | e' a f)) s (Y) = J\f(Y \ f) + Tr s+1 ^(-Mntfr s + 1)) (/ s+ i(Y, a + 1) + (50) 



P:(y,s + l) = X 1 |J^2, 
!£li;s+lel 2 



where we use notations as above for formula (l24"j) . and transformations similar to equalities (T5j 
and (l31~j) have been applied. 

Indeed, to set down a generator of mapping (I49p in componentwise form we observe that 
according to definitions (HTT) and ([3]), the following equality holds 

CO _ 

(e a AT(- | e-/)).(y) = J2 ~ Tr s+ i,...,s + n ( - ^(^/UW + 

n=0 

+ E E E(-- A/ "^(i 1 ,i 2 ))(e-7)|x 1 |(X 1 )(e-7)|x 2 |(X 2 )). 
Then in view of formulas (1411) and (H3l) we have 



(e n /)|Xi|(^l)(e a f)\X 2 \(X 2 ) = ^ 7j TTs+n+l,...,s+n+fc E I \7l fc \; /|^i|+fc-fci (^1; 

fc=0 ' fci=0 

s + n + l,...,s + n + k — k 1 )f\ X2 \+k 1 (X 2 , s + n + k — h + 1, . . . , s + n + k), 
and result we obtain 

(e°jV(- | e- a f)) s (l,...,s) = 

v ^ n | 



n\ ' fc!(n — A;)! 

n=0 fc=0 v y 



+ E 



k\ 



k x \(k-k x )\ 

P:(i,...,s+n-k)=x 1 \jx 2 hex 1 i 2 ex 2 fei=0 1V ±y 



X 



x/|Xi|+fe-fei(^i, s + n-fc+l,...,s + n- fci)/|x 2 |+ fcl (X 2 , s + n - k x + 1, . . . , s + n)) . 
Therefore the first term of this series is generator ([3]) of the von Neumann hierarchy 

h = (-K)fs(Y) + E £(-M*(n,ia))W*i)W*2). 

P:Y=X 1 (JX 2 heXx i 2 £X 2 
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and, as stated above, it coincides with the generator of the first order cumulant ( 14"5"1) . For the 
second term of series fl50l) we have 



I 2 = Tr s+1 ((-Af s+1 (Y, s + - (-K(Y))f s+1 + 

+ E E £ (-^*(m 2 )) w*o w**) - 
- E E E(- a/ '^(^i,^))(/ix 1 i + i(^i, s +i)/ix 2 i(x 2 ) + 

p-.Y=x 1 [jx 2 hex! i 2 ex 2 
+f\x 1 \(X 1 )f lX2l+1 (X 2 ,s + l))). 

Taking into account equality (TSSJ) in case of the set (Y, s + 1), it holds 

E E E(- A/ '^^'^))/ix 1 |(x 1 )/| X2 |(x 2 ) = 

P: (y,s+i)=Xi (J x 2 ueXi i 2 ex 2 

E EE (--Mnt(^l^ 2 ))/|X 1 |(^l)/|X 2 |(X 2 ) + 

P: (y,s+i)=X! u x 2 neii n y *2ex 2 f| Y 

+ E E (-MUn,s + l))f\ Xl \{X 1 )f\ X2 \{X 2 ) = 

p : (y,s + i) = x 1 \jx 2 , iieXipy 
s + iex 2 

= E E E(" A/ '^^'^))(/|y 1 | + i(F 1 , S + l)/| y2 |(r 2 ) + 
p : y=y 1 yy 2 ney « 2 ey 2 

+/|y 1 |(^i)/|y 2 | + i(^2, s + 1)) + J^HM 4 ' * + !)) E W^OfoiM, 

«ey p : (y, s + i) = XiU^2, 

i£Xi;s + leX 2 

where ^ p (y,s + i) = x 1 \jx 2 , is the sum over all possible partitions of the set (Y, s + 1) into 

s + 1 e x 2 

two mutually disjoint subsets X 1 and X 2 such that (s + l)£/i particle index belongs to set X 2 . 
As a result we obtain 

J 2 = Trs+i^i-Mrtii, s + l))/ s+1 + 

igy 

+ X)(-^U(<,« + 1)) £ feiW/pr,!^)), 

*ey p : (y,s + i) = XiU^2, 

» e Xi ; s + i e x 2 

i.e. this term coincides with the generator of second-order cumulant (1451) . 

In case of a two-body interaction potential other terms of series (150]) are identically equal to 
zero. This statement is a consequence of the structure of expansion fl50|) and of the fact that 
its third term equals zero. Indeed as a result of regrouping terms in the expression of the third 
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term we obtain 

I a = iTr s+M+2 ((-A/; +2 (y, 8 + 1,b + 2)- 2!(-JV <+ i(y, 8 + 1)) + (-W s (F)))/ s+2 + 
+ E EE (-Mnt^i,^))/!^!^!)/^!^)- 

P: (Y,s+l,s+2)=X 1 {JX 2 heXi f]y i2€X 2 f\Y 

" 2! E EE (-^(n^2))(f\x ll+ i(x u s + i)f lX2l (x 2 ) + 

p : (Y, s +i)=Xi u x 2 hex! n y i 2 ex 2 f] Y 
+f\x 1 \(X l )f l x 2 \+i(X 2 ,s + l)) + 

+ E EE (-Mint(h,i2W\x l \ +2 (X l ,s + l,s + 2)f lX2l (X 2 ) + 

p : k=xi u ^2 »i exi n y ^ax 2 n v 

+/|x 1 |(X 1 )/| X2 | +2 (X 2 , s + 1, 8 + 2) + 2\f ]Xll+1 (X 1 , s + l)f ]X2l+1 (X 2 , s + 2))) = 0. 

Thus, we conclude the validity of formula (1501) in case of a two-body interaction potential 
which describes the structure of the infinitesimal generator of mapping (j4*9l) in the general case. 

4.4 An existence theorem 

For an abstract initial- value problem of hierarchy (124p in the space the following theorem 

is true. 

Theorem 1. J/max„>i \ < (2e 3 ) _1 ; t/ien m case of bounded interaction potentials 

for t£l« solution of the Cauchy problem of the nonlinear quantum BBGKY hierarchy (24\)- 
( PSP determined by expansion (f^ ). If G n (0) G S^iT-L^) C fi 1 ^,,), ^ a strong (classical) 
solution and for arbitrary initial data G n (0) e £ 1 ('H n ) zs a wea/c (generalized) solution. 

Proof. It will be recalled that according to estimate (|48p . series (|44p converges under the con- 
dition: max n >i llGnfO)!! lf . < (2e 3 ) -1 . To prove that a strong solution of the nonlinear 

BBGKY hierarchy (|2"4"1) is given by expansion (jHj) we first differentiate it over time variable 
in the sense of a pointwise convergence on the space (H n ), i-e. for every function from the 
domain ^ 6 T>(H S ) C % s . Taking into account the group property of mapping (T4"9l) generated 
by expansion (j4"4"|) and properties of reduced nonlinear cumulants (|4"5l) of groups of operators 
©, for G n (0) E £l{U n ) C £ x (^ n ), n > 1, we obtain 

1 °° 1 

lim — ( V - Tr s+li ... jS+n (C/ 1+n (At; {y}, X \ F | G(t)) - G,(t, Y))^ s = (51) 
At->o Ac z — ' n! 

n=0 

= AT(y | G(t))ip s + Tr s+1 ^(-AT int (z, s + 1)) (G s+1 (t, Y,s + 1) + 

ieY 

P:(Y,s + l) = X 1 [JX 2 , 
i e X x ; s + 1 e x 2 

where expansion (j44|) is denoted by the symbol G s (t, Y). Since G n (0) G £\{(H n ) C -C 1 ('H„,), n > 
1, then using equality (15T|) . in the sense of the norm convergence in 2}(T-L n ) we finally establish 
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the validity of the equality 

1 °° 1 

lim Tr w | — ( V - Ti s+1 _ s+n U 1+n (t + At; {Y},X\ Y | G(0)) - 

n=0 

00 J 

- E ^ ^+1,..,*+^!+^; * \ y I G(o))) - 

n=0 

- (jV(y I G(t)) + Tr s+1 ^(-jVUCi, a + 1)) (G s+1 (t, V, a + 1) + 
+ Yl G lXll {t,X 1 )G lX2l {t,X 2 ))}\=0, 

P:(Y,s + l) = X 1 {JX 2 , 
i e X\ ; s + 1 e x 2 

which means that a strong solution of the nonlinear BBGKY hierarchy (|24l) is given by expan- 
sion (jUJ) in case of initial data from the subspaces £j\{^-L n ) C 2}{H n ), n > 1. 

Let us give a sketch of the prove that in case of arbitrary initial data G n (0) G £ 1 ('H n ), n > 1, 
expansion (jUJ) is a weak solution of the initial- value problem (I24p - (l2"5"|) . To this end we introduce 
the functional 

00 1 

(/, G(t)) = J2~ 'in y), (52) 

where / = (0, fx, . . . , f n , . . .) G fi (^ : W) is the finite sequence of degenerate bounded operators 
with infinitely times differentiable kernels with compact supports. For G n (0) G 2}{1-L n ) and 
f n G £o(% n ) functional ( 15"2"|) exists. 

We transform functional ( 152]) to the following form 

(/, G(t)) = (/, e a G(t | e- a G(0))) = (e a+ /, 0(t | e- Q G(0))), (53) 

where the operator o is defined by (l4"Tj) and on / s G £oCH s ) the operator a + is defined by the 
formula (an analog of the creation operator) 

(a+fUY) = Y,fs-i(Y\U)). 
i=i 

To differentiate obtained functional ( 153]) with respect to the time variable we use the corre- 
sponding result [12] of the differentiation of group ([7]) of the von Neumann hierarchy (JT|). As a 
result we derive that 

jSfMt)) -^^■■A^{Y)(e a+ f) s (Y)g(t,Y | e -G(0)) + 
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Taking into account the structure of expansion (jHJ) of the nonlinear quantum BBGKY hierarchy 
solution, for f s G £,q(% s ), s > 1, the following equality holds 

j t (f,G(t)) =£^ W (WW + E Mint(i,j)fs-i(Y\(j)))G s (t,Y) + 

s=0 ' i,j = l 

+ E EE E M nt (z 1 ,^ 2 )/ s (y)G| Xl |(t,X 1 )G|x 2 |(t,X 2 ) + 

P:Y=Xl U^2 iieXi ii€Xi i 2 ex 2 
s 

+ E ^U^U-itnO')) E G lXll (t,X 1 )G lX2l (t,X 2 ))y 

i,j = l P:Y = X 1 (JX 2 

i + j i £ Xi;j e x 2 

This equation means that in case of arbitrary initial data G n (0) e £ 1 ('H n ), n > 1, a weak 
solution of the initial- value problem (I24p - (l2"5"j) is given by expansion (JHJ). 

□ 



4.5 Remark: the nonlinear Vlasov hierarchy 

We give comments on the mean field asymptotic behavior [23] of constructed solution 
Let us suppose the existence of the mean field limit of initial state in the following sense 

lim||e s G^0)-£ n (0)|L._. = 0, n > 1. (54) 
Then there exists the mean field limit g s (t, 1, . . . , s), s > 1, of marginal correlation operators 



\mJe s G s (t)-g s (t)\\ 2HHg) = a, 8 >1, 
which is governed by the nonlinear Vlasov quantum hierarchy 

jg s (t, Y) = J2(-M(i))g s {t, Y) + Tr s+1 Y,{~M»t(h s + l))x (55) 



X 



(g s+1 {t,Y,s + l)+ Yl 9\x 1 \(t,X 1 )g lX2l (t,X 2 )), s > 1, 



P : (Y,s + l) = X 1 \JX 2t 

i e x 1 ■ s + 1 e x 2 



where notations similar to hierarhy (1241) are used. 

If initial data satisfies chaos property, then we establish 



lim ||e s C s (i)|| £l(Ks) = 0, s>2, (56) 



since solution expansions fl37|) for marginal correlation operators are defined by the (s + n)th- 
order cumulants as contrasted to solution expansions (122]) for marginal density operators de- 
fined by the (1 + n)£/i-order cumulants and in the consequence of the following formula on an 
asymptotic perturbation of cumulants of groups of operators [24] 



lin k\\\^s+n(t,l, ■ ■ ■ , s + n)f s+n \\ 1( = 0. 
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In case of s = 1 provided that ( 1541) we have 

lim \\eG x {t) - #i(£)|| £ i(k) = ° 3 



where for finite time interval the limit one-particle marginal correlation operator gi(t, 1) is given 
by the norm convergent on the space 2}(Ji) series 

<7i(M) = (57) 

~, * *n— 1 n 

= X) / d *i ■ - ■ / dt n ^ 2 ,„,^i(-t + ti,l)(-jV^(l,2))n^ 1 (-ti + t 2 ,ji)... 

"=0 £ £ ii=l 

n n n+1 n+1 

n "I - ^n; ^n 

which obviously coincides with iteration series of the Vlasov quantum kinetic equation [10]. For 
bounded interaction potential jT]) series (1571) is norm convergent on the space 2 1 (H) under the 
condition: t < t = (2 ||$|U(w 2 )|bi(0)|| £ i( W) ) . 

In view of the validity of limit (|56|) from the Vlasov nonlinear quantum hierarchy (|55j) we also 
conclude that limit one-particle marginal correlation operator (1571) is governed by the Cauchy 
problem of the Vlasov quantum kinetic equation 

^ gi (t, 1) = -Af(l)gx(t, 1) + Tr 2 (-AUl, 2))gi(t, l) gi (t, 2), (58) 
at 

and consequently for pure states we derive the Hartree equation. 

Thus, the nonlinear Vlasov quantum hierarchy (|55|) describes the evolution of initial corre- 
lations. 



5 Conclusion 

In the paper the origin of the microscopic description of non-equilibrium correlations of quantum 
many-particle systems obeying the Maxwell-Boltzmann statistics has been considered. The 
nonlinear quantum BBGKY hierarchy ([21]) for marginal correlation operators was introduced. 
It gives an alternative approach to the description of the state evolution of quantum infinite- 
particle systems in comparison with quantum BBGKY hierarchy for marginal density operators 
[THIHI]. The evolution of both finitely and infinitely many quantum particles is described by 
initial-value problem of the nonlinear quantum BBGKY hierarchy (T2^1) and in case of finitely 
many particles the nonlinear quantum BBGKY hierarchy is equivalent to the von Neumann 
hierarchy <^). 

A nonperturbative solution of the nonlinear quantum BBGKY hierarchy is constructed in 
the form of expansion (|44|) over particle clusters which evolution is governed by corresponding- 
order cumulant (|4"5I) of the nonlinear groups of operators generated by solution flTJ) of the von 
Neumann hierarchy (Q. We established that in case of absence of correlations at initial time the 
correlations generated by the dynamics of quantum many-particle systems (1371) are completely 
determined by cumulants (|Hj) of groups of operators dHj). 
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Thus, the cumulant structure of solution (J7J) of the von Neumann hierarchy (TjQ) induces the 
cumulant structure of solution expansion (I44p of initial- value problem of the nonlinear quantum 
BBGKY hierarchy ([M]). 

We emphasize that intensional Banach spaces for the description of states of infinite-particle 
systems, which are suitable for the description of the kinetic evolution or equilibrium states, are 
different from the exploit spaces [H], |19j . Therefore marginal correlation operators from the 
space of trace-class operators describe finitely many quantum particles. In order to describe 
the evolution of infinitely many particles we have to construct solutions for initial data from 
more general Banach spaces than the space of sequences of trace-class operators. For example, 
it can be the space of sequences of bounded translation invariant operators which contains the 
marginal density operators of equilibrium states [25J. In that case every term of the solution 
expansion of the nonlinear quantum BBGKY hierarchy (J33J) contains the divergent traces, which 
can be renormalized due to the cumulant structure of solution expansion (T4"5"l) . 

The mean field asymptotic behavior of constructed solution (jHJ) is governed by the nonlinear 
Vlasov quantum hierarchy (1551) . In such approximation this hierarchy describes the evolution 
of initial correlations and in case of its absence the nonlinear Vlasov hierarchy (|55|) is equivalent 
to the Vlasov quantum kinetic equation (|5"81 . 

Following to the paper [12] the obtained results can be also generalized on many-particle 
systems obeying the Fermi-Dirac and Bose-Einstein statistics (1321 . 
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